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O ! Abstract 
(N 

^jq' We study the finite-size effects on tlie deconfinement pfiase transition (DPT) of liot and / or dense 
^3 ■ 

■ hadronic matter, using a simple tliermodynamic model based on the assumption of coexistence 

. of confined and deconfined phases in a finite volume, with the two-phases matter Equations of 



p. 

> 

X 



requirement. A problem arises in the limit of small QGP volumes when using the usual color-singlet 
partition function (CSPF) derived in the saddle point approximation. To avoid this problem, 



; State (EoS). For the QGP, we consider a partition function (PF) with the exact color-singletness 

m 

i we have then proposed a method for calculating a suitable CSPF which allows us to accurately 

o ■ 

calculate physical quantities describing well the DPT at finite volumes like the order parameter, 
the energy density and the entropy density. We show that in the limit of infinite volume, these 
D , thermodynamic quantities exhibit a discontinuity at a critical temperature Tc(oo) if the transition 

is temperature driven, or chemical potential /Xc(oo) if it is density driven. In a finite size system, 
^ ■ all singularities are smoothed out over a broadened critical region, shifted from the critical point 

position in the thermodynamic limit. At the level of the first derivatives of these thermodynamic 
quantities, the specific heat and the susceptibility show delta function singularities at the critical 
point in the thermodynamic limit. In a finite volume, these delta singularities are smeared into 
finite peaks of widths 5T(y) or 5fi{V), with the maxima of the peaks occuring at pseudo-critical 
points Tc{V) or HciV). For the temperature driven DPT, an analysis of the finite size scaling 
behavior at criticality of these maxima as well as of the width of the transition region and the 
shift of the transition temperature allows us to determine the critical exponents characterizing the 
deconfinement phase transition. The obtained results are in good agreement with those predicted 
by other studies for a first-order phase transition. 
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I. INTRODUCTION 



Phase transitions in statistical physics are known to be infinitely sharp only in the ther- 
modynamic limit {V — > oo). Only in this limit is the thermodynamical potential or any 
of its derivatives singular at the critical point. However, real systems and the systems we 
simulate are finite. For instance, the reaction zone in which the formation of a quark-gluon 
plasma (QGP) in a relativistic heavy ion collision is assumed to take place is limited, with a 
width of about 1 fermi, and a transverse area of a radius comparable to the radii of the col- 
liding nuclei, i. e., V ^ A^^^fm^. Also, lattice QCD calculations are necessarily performed 
in finite lattices. 

In general, finite size effects lead to a mixed phases system and a rounding of the tran- 
sition. The main care is then how to sign a possible phase transition in a finite system. 
It turns out that we can extract the true critical behavior of infinite systems from calcu- 
lations on finite systems, by studying how some characteristic thermodynamical quantities 
vary with the size of the system, namely by a Finite Size Scaling (FSS) analysis. Studies in 
statistical physics have shown that despite the apparent diversity in the underlying struc- 
ture of systems undergoing phase transitions, these latter take place with some universal 
global behavior, depending only on the range of interaction of the forces at play and the 
dimensionality of the problem. This means the singular behavior of some characteristic 
observables near criticality is identical for many systems when appropriately scaled. This 
universal behavior of finite systems at criticality is characterized by indices called critical 
exponents. The determination of critical exponents has long been one of the main interests 
for both analytical calculations and numerical simulations. 

In the present work, we study the finite size effects on the deconfinement phase transition 
(DPT) when temperature driven and when density driven, within a simple thermodynamical 
model used in [0], based on the assumption of the coexistence of confined and deconfined 
phases in a finite volume, and using the equations of state (EoS) of the two phases. For the 
QGP, we consider a partition function with the exact color-singletness requirement. This 
latter is usually derived using the saddle point approximation, and gaussian approximation 
which break down at VT^ << 1 [0, |^, ^ . To avoid the problem at small temperature and / 
or volumes, we have then proposed a method for calculating a suitable CSPF, expanded in 
a power series of the QGP volume fraction, which allows us to accurately calculate physical 
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quantities describing well the DPT at finite volumes like the order parameter, the energy 
density and the entropy density. We show that in the limit of infinite volume, these thermo- 
dynamic quantities exhibit a discontinuity at a critical temperature Tc(oo) for a temperature 
driven DPT. In a finite size system, no singularity occurs and the variations of these quan- 
tities are perfectly smooth over a broadened critical region, shifted from the critical point 
position in the thermodynamic limit. 

At the level of the first derivatives of these thermodynamic quantities, the specific heat 
and the susceptibility show delta function singularities at the critical temperature in the 
thermodynamic limit. In finite volumes, these delta singularities are smeared into finite 
peaks of widths 6T{V), with the maxima of the peaks occuring at pseudo-critical temper- 
atures Tc{V). An analysis of the finite size scaling behavior at criticality of these maxima 
as well as of the width of the transition region and the shift of the effective transition tem- 
perature Tc(y) relative to the true one Tc(oo) allows us to determine the critical exponents 
characterizing the deconfinement phase transition. 



II. COLOR-SINGLET PARTITION FUNCTION OF THE QGP 

The partition function for a color-singlet quark-gluon plasma contained in a volume Vqgp, 
at temperature T and quark chemical potential /i, is determined by |0, ||, ^: 

Z(T, Vqgp, /x) = ^e^^ £^ £^ d (|) d (^^^ M(^, V^)Z(T, Vqgp, /i; ^, (1) 

where the factor exp(— '^^^°^ ) accounts for the vacuum pressure B, M{(p,ijj) is the weight 
function (Haar measure) given by: 

M(^,^)= (smQ(^ + |))sin(|)sinQ(^-|)^^ , (2) 

and Z the generating function defined by: 

Z(T, Vqgp, fi; ^, ^) = Tr [exp (^o - ft (^Ng - iV^)) + i^h + ^V'^'s)] (3) 

where (3 = — with the units chosen as: ks = h = c = 1, Hq \s the free quark-glue 



Hamiltonian, Nq \Nqj denotes the (anti-) quark number operator, and /a and are the 
color "isospin" and "hypercharge" operators respectively. 
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The generating function Z (T, Vqgp, IJ>; (p, ip) can be factorized into the quark contribution 
and the glue contribution as: 

Z {T, Vqgp, ip, ^) = Zquark {T, Vqgp, ip, ip) Zgi^e {T, Vqgp; ip, ^) , (4) 
where the quark contribution is given by: 



Zquark {T, Vqgp, H', V') = ©Xp 



-T'VQGpdQ J2 [ ^ 

q=r,g,b 



TT 



+ 2V . 



(5) 



with q = r, b, g the color indices, (Iq = 2Nf counts the spin-isospin degeneracy of quarks 
and the angles aq being determined by the eigenvalues of the color charge operators in eq. 
(3): 



"^=2 + 3' "^ = "2 + 3' "^ = -T' 



(6) 



and the glue contribution is given by: 
Zgiue {T, Vqgp; V') = exp 



.9=1 \ ^ 



a„ — n\ 1 I ag — Tx 



2 V TT 



with (Ig — 2 the degeneracy factor of gluons and ag {g — 1, ..A) being: 

ai = ttr — ag, a2 = ag — ab, a^ = ab — ar, a^ = 0. 
Let us write the generating function as: 

Z {T, Vqgp, /x; cp, ^) = e^^'^-^'^^^'^'^), 

with: 

/ , /^N ^%21 , 16,, n^dn V- f . A' 

g((p, ^, -) = —i—dQ + —dc) + > < -1 + — - 1 

^^^'^'T^ 12^30 ^ 15 ^ 12 2 ^ I I 7r2 i 

q=r,b,g \ / 

2 J 4 



(7) 



(8) 



(9) 



12 2 



E 

9=1 



16 J , 




' 12 2 


("a - ^] 
7r2 





(10) 



then eq. (1) can be written on the form: 



ZiT,f,,VQGp) = Ae-^(^^^o{^)-f) p p di^^/2)d{^/?,)M{^,^) 



M^,^,m-go{^))xVT^ 



(11) 
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where go (^) is the maximum of g{ip, tp, ^) for (p,ip & [— tt, +7r], given as: 

vr^ 21 16 tt^c/q 3yU^ 3yU^ , . 

ItJ = 12^30^'^ + 15^^^ + -IT^^^ + 2^^' ^^^^ 

and X the QGP volume fraction defined by: x = 1 — f), since assuming the phases coexistence 
model, the hadronic gas and QGP phases occupy fractional volumes Vhg = and Vqgp = 
(1 — f)) V, characterized by the parameter f) lying in the range — 1, where f) = 1 corresponds 
to a pure hadron phase and f) = to a pure QGP phase. 

When VT'^ ^ 1, we can evaluate the integral in eq. (11) by the saddle point approxima- 
tion around the maximum of the integrand at (0, 0). After some calculation, eq. (11) can be 
put on the form [^: 

Z{T,^^Vx)^^ ' / / rf^rf^M(°'°)(^,V.)e-'"'"'(-^^^^^'^), (13) 

with: a (f ) = + ^ i^^) ' expansions of M and {g - ^o) 

respectively, around (v^,'?/') = (0,0), given by: 



127r2(a(^))VT3 V8 9 3 
Let us write the exponential in the integral of eq. (13) for VT^ ^ 1 as: 

V 7r(a(^))'/'v^ V4 

+ ^ + (15) 



127r2(a(^))VT3 V8 9 3 

and do the expansion: 



5^ AZlA^l_p^^. ; (16) 
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then we obtain for the CSPF of the QGP: 
WT,.,Vx)^^ ^g"^-^^-12.^(a(.))VT3g^^-^" ), (17) 
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such as the coefficients aj and (3j of the two series are given by: 



It can easily be shown that the series in eq. (17) are convergent, and can be truncated 
at a rank j above which the contribution of the terms ajx^ to the sums is negligible. 

III. FINITE SIZE EFFECTS 

To study the effects of volume finiteness on the DPT within the phases coexistence model 
(PCM), let us examine the behavior of some thermodynamic quantities in the vicinity of 
the critical point. The quantities of interest are the order parameter, which is simply in this 
case the hadronic volume fraction, the energy density and the entropy density. The mean 
value of any intensive thermodynamic quantity of the system can be calculated by: 

/ ((1 - ()) ((1 - l))V) + () Ah iW)) Z (()) di) 

{A{T, V)) = ° , (19) 

IZ{i)) di) 



where Aqgp (f)) and Ahg (f)) are contributions relative to the individual QGP and hadronic 
gas phases, respectively, and Z {i)) the total PF of the mixed system, which factorizes into 
the QGP and hadronic gas partition functions as: Z (i)) — Zqgp (1)) Zhg (f)) ■ 

For the hadronic gas phase, let us consider the partition function of a pion gas given by: 

Zhg = e-w^^^- . (20) 
The mean value of the order parameter is then given by: 

/dxexp(((g + ^o(^))T3-f)xl^) (Er=o«.^^+^- ' 

<t){T,pi,V)>^l 



127r2 {a{§)) VT\ 



1 / 7Y°° d-x^^^ 

/d,xexp(((f + ^o(^))T3-f)xV^) Er=o«.^^- ' 







127r2 {a{§)) VT\ 
(21) 
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and those of the energy and entropy densities can be written, respectively, as: 

< s{T, V) >= CHG + eo + (eQGP - chg) < x{T, V) > (22) 

< s{T, V) >= shg + Sq + (sqgp - shg) < x{T, n,V) >, (23) 

with: 

< x{T,fi,V) >= 1- <i){T,fi,V) > 

(24) 

As it is known, the deconfined quark matter state can be obtained at extreme conditions 
of temperature and/or density, i. e., either by raising temperature or chemical potential. 
The deconfinement phase transition is then temperature driven or density driven. 

We'll examine in the following the FSE for a temperature driven DPT, at a van- 
ishing chemical potential (/x = 0), and for a density driven DPT at a fixed tempera- 
ture (T = lOOMeV) , considering the two lightest quarks u and d {Nf = 2), and using 
^1/4 — i4^^MeV for the bag constant. 

Fig.(l-a) illustrates the variations of the order parameter, the energy density normalized 
by T^, and the entropy density normalized by T^, versus temperature for different system 
sizes. The curves show that in the limit of an infinite volume, the three quantities exhibit 
a sharp discontinuity at a transition temperature Tc (cx)), reflecting the first order character 
of the transition. For small size systems, the transition is smoothed out over a range of 
temperature ST{V), around a pseudo-critical temperature Tc{V), shifted from the true 
transition temperature Tc (oo). The broadening of the transition region, as well as the shift 
of the critical temperature get larger, smaller is the volume. 

In the same way as for temperature, we can deal with the density driven phase transition, 
at a fixed temperature (T = lOOMeT^). The variations of the order parameter, the energy 
density and the entropy density versus chemical potential are presented for various system 
sizes on Fig. (1-b). The first-order character of the transition can, also in this case, clearly 
be seen at the large volumes limit. In small systems, the transition is perfectly smooth over 
a region of chemical potential of width S/j, (V), and the effective critical chemical potential 
/ic (V) is shifted away from the true one //c (oo) . 
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IV. FINITE-SIZE SCALING ANALYSIS 



A. Finite size scaling and critical exponents 

In the thermodynamic hmit, phase transitions are characterized by the appearance of 
singularities in some second derivatives of the thermodynamic potential, such as the sus- 
ceptibility X cLnd the specific heat c. For a first order phase transition, the divergences are 
typically 5-function singularities, corresponding to the discontinuities in the first derivatives 
of the thermodynamic potential. 

For temperature driven phase transitions, in finite volumes, the singularities in x {T, V) 
and c (T, V) are rounded over a range of temperature 6T (V) . The peaks occuring at a 
pseudo-critical temperature Tc {V) may be shifted away from the true critical temperature 
Tc(oo). The width of the transition region, the shift of the critical temperature and the 
maxima of the peaks of the specific heat and the susceptibility, represented in Fig. (2), show 
a scaling behavior at criticality, characterized by critical exponents as: 

' 5T (V) ~ V-^^ 

TTiv) = niv)-ni^)r^v-'- 

< [2b) 

c^^^ {V) ^ V""^ 

. xT''{y)-y^^- 

The finite size scaling (FSS) analysis is used to recover the critical exponents, and it has 
been shown that for a first order phase transition , the finite size quantities 5T (V), TxiV) 
scale as: V^^ and the maxima c^^^ iV) and Xt^^ (^) scale as: V . The critical exponents 
Oti At, otT and 7t are then all equal to the dimensionality. 



B. Numerical determination of the critical exponents 

1. The susceptibility critical exponent 

For the temperature driven phase transition, the susceptibility is defined to be the first 
derivative of the order parameter with respect to temperature, i. e., 

„r,n^^^*§^. (26) 
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It is clear in this case from Fig. (3-left) that the delta function singularity of the susceptibihty 
in the thermodynamic limit is smeared, in a finite volume, into a finite peak of width 5T {V), 
with the maximum of the peak | xt (Y) occuring at the pseudo-critical temperature 
Tf.{y). The plot of the maxima | xt |™^^ iY) versus volume is illustrated in Fig. (3-right), 
and the linearity of the data with V can clearly be noted. A numerical parametrization 
with the power- law form: | xt |"^^ iV) ~ V'^'^ , gives the value of the susceptibility critical 
exponent 77- The determination of the location of the maxima being done using a numerical 
method with an error on the temperature and on the maxima | xt iy)-i this yields 
a systematic error on the exponent 7t. The systematic error is given for all the critical 
exponents determined in the following. The result for the susceptibility critical exponent 
and the estimated systematic error are then: 7^ = 0.99 ± 0.04. 

2. The specific heat critical exponent 

The specific heat density representing the smeared delta function of the latent heat is 
defined as: 

e being the energy density of the system. The variations of the specific heat density with 
temperature are presented for different sizes in Fig. (4-left) , which shows the rounding of the 
delta function singularity of c (T, V) in finite systems into a finite peak of width 5T (V) and 
height c™^^ iy ) . For decreasing volume, the width get larger while the height of the peak 
decreases. The data of the maxima of the specific heat c™^^ {V) are fitted to the power-law 
form: c™*^^ {V) ~ V^"^ in Fig. (4- right), and the obtained specific heat critical exponent is: 
ar = 0.99 ± 0.04. 

3. The shift critical exponent 

For the study of the shift of the transition temperature rr(T^) = Tc (V) — Tc (00), we need 
to have the effective transition temperature in a finite volume. It is usually defined to be 
the temperature at which the rounded peaks of the susceptibility and the specific heat reach 
their maxima. 



9 



Fig. (5) illustrates the results for the shift of the critical temperature plotted versus 
inversed volume, and shows the linear character of the variations. The shift critical exponent 
obtained from a fit to the form: rriV) ~ F"^^, is: At = 1.0085 ± 0.0009. 



4- The smearing critical exponent 

The width of the transition region can be defined by the gap: ST{V) — T2{V)—Ti{V) with 
Ti(V) and T2{V) the temperatures at which the second derivative of the order parameter 
reaches its maxima, or in other terms the temperatures at which the third derivative of the 
order parameter vanishes, i. e., 

Q3 



{h{T,V)) 



= 0. (28) 

ri(v),T2(y) 



Fig. (6-left) illustrates the variations of the second derivative of the order parameter 
with temperature for various sizes, and shows that the gap between the two extrema, which 
represents the broadening of the transition region, decreases with increasing volume. The 
results for the widths 6T {V), plotted in Fig. (6-right) vs volume, were fitted to the power law 
form: ST {V) ~ V~^'^, and the obtained smearing critical exponent is: 9t = 1.016 ± 0.007. 



V. CONCLUSION 



Our work has shown the influence of the finiteness of the system size on the behavior 
of thermodynamical quantities near criticality. The sharp transition observed in the ther- 
modynamic limit, signaled by discontinuities in the first derivatives of the thermodynamic 
potential at a critical temperature or chemical potential, is rounded off in finite volume, and 
the variations of the thermodynamic quantities are perfectly smooth on the hole range of 
temperature or chemical potential. For the temperature driven DPT, a finite-size scaling 
analysis of the behavior of the width of the transition region ST (V), the shift of the pseudo- 
critical temperature relative to the true one tt (V) — Tc {V)—Tc (oo), and the maxima of the 
rounded peaks of the susceptibility Xt^"^ (^) and the specific heat c^'^"^ (V) near criticality. 
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shows their power-law variations with the volume: 

5T {V) ~ V'^^T 

TTiV) = T, iV) - T, (oo) ~ V-^^ 
c^^^ [V) ~ 

characterized by the scaling critical exponents 6t, At, and jt- Numerical results for 
these critical exponents have been obtained, and the associated systematic errors, resulting 
from the numerical method of calculating the maxima Xt**^ (^) ^^"^ '^t^^ (^) ^^"^ their 
localizations as well as the gap 6T(y) = T2(y) — TiiV), have been estimated. Our results 
for the critical exponents are in good agreement with the analytical ones; O'p — Xx — ttT — 
7t = 1 found in a parallel work [^. These results are characteristic of the first order phase 
transition, as predicted by the FSS theory. 

A further FSS analysis is viewed using an exact color-singlet partition function of the 
QGP, derived without any approximation, with an appropriate numerical treatment. Also, 
the density driven DPT can be studied, and similar critical exponents associated to the shift 
of the critical chemical potential [fic (V) — fic (oo)], the width of the transition region 6fi {V) 
and the maxima of the rounded peaks of the susceptibility x™'^^ (^) specific heat 

c^^^ {V) (relative to the derivatives with respect to fi) can be determined. 
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Figure 1: Plots of (a) order parameter (bottom), energy density normalized by (middle) and 
entropy density normalized by (top) versus temperature at = 0, and of (b) order parameter 
(bottom), energy density (middle) and entropy density (top) versus chemical potential at T = lOOMeV, 
for different system sizes. 
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Figure 2: Illustration of the critical behavior of the susceptibility x {T, V), the specific heat c (T, V) 
and the second derivative of the order parameter dx/dT. 




Figure 3: (Left) Susceptibility x {T, V) as a function of temperature for varying volume, and (Right) 
linear fit of the results for the maxima of the susceptibility vs volume. 




Figure 4: Variations of (left) the specific heat density with temperature for different system volumes 
and of (right) the maxima of the specific heat density with the volume. 




Figure 5: Plot of the shift of the critical temperature vs inversed volume. 




Figure 6: (Left) Illustration of the second derivative of the order parameter ^ ^g^'^^^ vs 
temperature, taken at different volumes. (Right) Data of the width of the temperature region over 
which the transition is smeared, fitted to a power-law of the inversed volume. 



